Spin-orbital Texture in Topological Insulators 
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Relativistic spin-orbit coupling plays an essential role in the held of topological in- 
sulators and quantum spintronics. It gives rise to the topological non-trivial band 
structure and enables electric manipulation of the spin degree of freedom. Because of 
the spin-orbit coupling, rich spin-orbital coupled textures can exist both in momentum 
and in real space. For three dimensional topological insulators in the Bi2Se3 family, 
topological surface states with p z orbitals have a left-handed spin texture for the upper 
Dirac cone and a right-handed spin texture for the lower Dirac cone. In this work, we 
predict a new form of the spin-orbital texture associated with the p^ and p y orbitals. 
For the upper Dirac cone, a left-handed (right-handed) spin texture is coupled to the 
"radial" ("tangential") orbital texture, whereas for the lower Dirac cone, the coupling 
of spin and orbital textures is the exact opposite. The "tangential" ("radial") orbital 
texture is dominant for the upper (lower) Dirac cone, leading to the right-handed spin 
texture for the in-plane orbitals of both the upper and lower Dirac cones. A spin- 
resovled and photon polarized angle-resolved photoemission spectroscopy experiment 
is proposed to observe this novel spin-orbital texture. 

PACS numbers: 71.20.-b,73.43.-f,73.20.-r 



I. INTRODUCTION 



Three-dimensional topological insulators (TIs) are new 
states of quantum matter with helical gapless surface 
states consisting of odd number of Dirac cones inside 
the bulk band gap protected by time-reversal symmetry 
(TRS). 1 - 4 The underlying physical origin of the topo- 
logical property of TIs is the strong spin-orbit coupling 
(SOC), which plays a similar role as the Lorentz force in 
the Quantum Hall state. Due to the SOC interaction, the 
spin and momentum are locked to each other, forming a 
spin texture in the momentum space for the surface states 
of TIs 5-7 . The spin texture has been directly observed 
in the spin-resolved angle-resolved photon emission spec- 
troscopy (spin-resolved ARPES) 8 ~ 12 . The spin texture 
gives rise to a non-trivial Berry phase for the topological 
surface states and suppresses the backscatterings under 
TRS, leading to possible device applications in spintron- 
ics. 

Besides the spin texture, it has also been shown re- 
cently that the atomic p orbitals of the Bi2Se3 family of 
topological insulators form a pattern in the momentum 
space, dubbed as the orbital texture, for the topological 
surface states. 13,14 In this work, we predict a coupled 
spin-orbital texture for the topological surface states. 
Based on both the effective k- p theory and ab-initio cal- 
culations, we find, besides the usual locking between the 
electron spin and the crystal momentum, the spin tex- 
ture is also locked to the atomic orbital texture, which is 
dubbed as "spin-orbital texture". We show that p z or- 
bitals have left-handed spin texture for the upper Dirac 
cone and right-handed spin texture for the lower Dirac 
cone, sharing the same feature as the total spin texture 
of the surface states. In contrast, the in-plane orbitals 





FIG. 1. (color online) a, b, The tangential orbital texture 
with the right-handed helical spin texture (a) and the radial 
orbital texture with the left-handed helical spin texture (b) 
for the upper Dirac cone. 



(p x and p y orbitals) reveal more intriguing features: for 
the upper Dirac cone of surface states, a "radial" orbital 
texture is coupled to a left-handed spin texture and a 
"tangential" orbital texture is coupled to a right-handed 
spin texture. For the lower Dirac cone, the coupling be- 
tween spin and orbital textures is exactly opposite. An 
electron spin-resolved and photon polarized ARPES ex- 
periment is proposed to observe this novel spin-orbital 
texture of the surface states of TIs. 



II. 



EFFECTIVE THEORY OF THE 
SPIN-ORBITAL TEXTURE 



The surface states of TIs are described by the Dirac 
type of effective Hamiltonian 5,15 



H smi {k x , k y ) = hv f (a x k y 



(1) 



with the Fermi velocity Vf and Pauli matrix a. The 
salient feature of this effective Hamiltonian is the "spin- 
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FIG. 2. (color online) a, The Dirac cone of Bi2Se3 on the surface with the normal direction [0001] with the spin texture marked 
by blue arrows, b, c, The projection of p z orbital and the related in-plane spin texture for upper (b) and lower (c) Dirac 
cones. More red means more p z character. The red arrows represent the in-plane spin texture related to the p z orbitals. The 
insets are the schematics of the spin texture marked by green arrows. 



momentum locking" , which means for a fixed momentum 
k, the "spin" , denoted by the Pauli matrix er, has a fixed 
direction for the eigenstate of the Hamiltonian. Since the 
"spin" is always perpendicular to the momentum, we can 
introduce a helicity operator, defined as h = ■ (k x a) 
which commutates with the Hamiltonian, to determine 
the handness of the "spin" texture. For the upper Dirac 
cone of surface states, the helical operator h = — 1, lead- 
ing to a left-handed "spin" texture in the momentum 
space while for the lower Dirac cone, h = 1 yields a a 
right-handed "spin" texture. However, one should note 
that here the "spin" is not the real spin, but the total 
angular momentum J = S + L, which is a combination 
of the real spin S and the orbital angular momentum L 
due to SOC. Consequently, the basis of the surface effec- 
tive Hamiltonian (1) are denoted as |^j z= ±|) with the 

lower indices ±i representing the total angular momen- 
tum along z direction. In order to understand what is 
the texture for the real spin S, it is necessary to write 
down the explicit form of the basis wavefunction l^^i). 

The form of the basis |\l/±i) can be constructed by 



symmetry considerations. Generally the basis l^^i) de- 
pends on the momentum k and we can expand it up to 
the first order in k as |* ± i) = I^YA) + Here we 

2 ^ ± 2 _ ±2 

are only interested in the p orbitals of Bi and Se atoms in 
the topological insulator Bi2Se3 and can decompose the 
zeroth-order wavefunction as 

l*±4> =EKa|a i ^ I ta)>+«0,a|«,P±,l(t)>] (2) 
a 

and the first-order wavefunction as 

l*±i) = 53 ± fc ±(«wi,a|a,p T ,t (I)) +ivi. a \a, Pz ,l (t)>) 

(3) 

=F iwi,a.k^\a,p±, J \ {{))] 

where k± = k x ± ik y , \ t) and | J,) denote the spin, \p z ) 
and \p±) = T^7| (\Px) ± APy)) denote different p orbitals, 
and a denotes indices other than the spin and orbital, 
such as atom indices. Here i«o(i),m v o(i),a an d w\^ a are 
material-dependent parameters. By comparing with the 
ab-initio calculations, we find that we can take them to 
be real. I^i) and are related to each other by 

TRS. The expressions of the basis (2) and (3) can be 
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FIG. 3. (color online) a, b, c, d, The p x projection on the states of upper (a) and lower (c) Dirac cones, and the p y projection 
of upper (b) and lower (d) Dirac cones. More red means more p x character in (a) and (c), and more red means more p y 
character in (b) and (d). The red arrows indicate the in-plane spin texture related to the orbitals. The insets are the schematics 
of the spin texture. 
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substituted into the eigen wavefunctions of the Hamil 
tonian (1), |$±) = -±= ±ie- l8 "\^i) + , yield: 

the following forms of the wavefunctions 

|$+) = ^2[{uo, a - v lyCc k)\a,p z ,te) 



1 

- -j={v Q , a - u 1>a k - w lja k)\a,PrA( 
+ -j=(v<o ta - u liCt k + Wx, a k)\a,p t ,U 



)(4) 



I*-) = 2j[(uo,a + vi >a k)\a,p z ,U) 



1 

+ -j=(v , a + Ui t0l k + Wi >a k)\a,p r ,U) (5) 

- -^={v , a + u lya k - w lya k)\a,p t ,te)}- 

Here | te (U)) = t^+R^-^ | t) + U)) stands for the 
left-handed (right-handed) helical spin texture and \p r ) = 
cos 9 k \p x ) +sm9 k \p y ), \p t ) = -sm9 k \p x ) + cos 9 k \p y } are 
the radial and tangential orbital textures, as shown in 



Fig. la,b, respectively. From the expressions (4) and (5), 
we can clearly see that \p z ) orbital is coupled to the left- 
handed spin texture | te) for the upper Dirac cone and 
the right-handed spin texture | 4-e) for the lower Dirac 
cone. Furthermore, for the upper Dirac cone, the radial 
orbital texture \p r ) is always coupled to the left-handed 
spin texture | \g) and the tangential orbital texture \pt) 
is always coupled to the right-handed spin texture | \.g). 
The situation is exactly opposite for the lower Dirac cone. 

The expressions (4) and (5) are the main analytical 
results of this paper, which show explicitly the spin- 
orbital texture. To confirm our analytical results, ab- 
initio method is adopted to calculate the projection of 
surface states on the spin and orbital basis, defined by 
the quantity 



D 



(*±|(|Pi)<Pi|«IS,)|*d 



(6) 



where pi — p x ,Py,Pz for the three p orbitals, sq = 12x2 
denotes the charge part and s x ^ y ^ z denote the three Pauli 
matrices for the spin. In the following, we will compare 
the analytical calculation of the intensity with ab- 
initio calculations. 




FIG. 4. (color online) a, b, The tangential orbital texture with the related in-plane spin texture for the upper Dirac cone (a) 
and the radial orbital texture for the lower (b) from ab-initio calculations. More red means more p x character, and more blue 
means more p y character. The red arrows represent the total in-plane spin texture related to p x and p y orbitals. c, The p x 
and p y projections onto upper and lower Dirac cones. The solid curves are the p x and p y projections for the upper Dirac cone 
at energy level O.lOeV, and the dashed for the lower Dirac cone at the energy — 0.07eV. The basic feature is the tt period which 
exactly agrees with the prediction by the effective model, d, Orbital polarization P Px . The energy of Dirac point is shifted to 
be zero. The positive value of the orbital polarization represents the radial orbital texture, and the negative value represents 
the tangential orbital texture. In order to plot more visually, the P Pa . for the upper Dirac cone is inversed, marked by the red 
color. P P:c is exact zero at Dirac point, which indicates the transition point between the tangential and radial orbital textures. 



III. SPIN-ORBITAL TEXTURE FROM 
AB-INITIO CALCULATIONS 

The Vienna Ab-initio Simulation Package (VASP) 16 ' 17 
is employed to carry out ab-initio calculations with the 
framework of the Perdew-Burke-Ernzerhof-type (PBE) 18 
generalized gradient approximation (GGA) of density 
functional theory 19 . Projector augmented wave (PAW) 
pseudo-potentials are used for all of the calculations in 
this work 20 . 10x10x10 and 10x10x1 are used for k-grid 
of bulk and free-standing calculations, respectively. The 
kinetic energy cutoff is fixed to 450eV. 6 quintuple layers 
(QLs) are fixed in the supercell for free-standing calcu- 
lations, and the thickness of vacuum is taken to be 50A. 
The lattice constant and the atomic parameters are di- 
rectly obtained from experiments. SOC is included with 
the non self-consistent calculation. In order to compare 
with the result of ARPES experiments, the projections 



of all the orbitals are only for the first Se and Bi atoms 
on the top surface of the free-standing model. 

The surface states of Bi2Se3 consist of a single Dirac 
cone at Y point on one surface inside the large bulk band 
gap (~ 0.3eV) 5,21 , which provides an ideal material to 
study the coupling of spin and orbital textures of sur- 
face states. As the starting point, we compare the bulk 
band structure of Bi 2 Se 3 with the previous calculation 5 
and find good agreements. The surface states are ob- 
tained from the calculation of a free-standing structure 
with the normal direction [0001], as shown in Fig. 2a. 
The blue arrows represent the spin texture, where the 
spin is mainly lying in plane near the Dirac point. The 
spin texture is left-handed for the upper Dirac cone and 
right-handed for the lower one, the same as the total 
angular momentum texture. To understand the under- 
lying physics, we calculate the spin texture for different 
atomic orbitals. For p z orbitals, a left-handed helical 
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spin texture is found for the upper Dirac cone and a 
right-handed texture for the lower Dirac cone, as shown 
in Fig. 2b, c. The schematic of the spin texture are 
shown in the inset. Here the background color indi- 
cates the projection of p z orbitals, which is isotropic, 
and the red arrows represent the corresponding in-plane 
spin texture. The spin texture of p z orbitals can be re- 
produced with the expressions [D^^^ , Dp x , ><T J = 
iEaKa T wi, Q fc) 2 [sin#fc, — cos#fc, 0] with '±' for the 
upper and lower Dirac cone and k = (k,9k) in the polar 
coordinate. 

The spin textures for in-plane orbitals are shown in 
Fig. 3a,b for the upper Dirac cone and in Fig. 3c, d for 
the lower Dirac cone, respectively. We find that the asso- 
ciated spins for p x and p y orbitals don't rotate clockwise 
or anti-clockwise around the Dirac point as in the case 
of p z orbitals, but instead, they take the form 

l D L*> D Ly> D P*J = T]T%[sin0 fe ,cos# fc ,O] (7) 



D ± D 1 



±^^[sin0 fe ,cos0 fe ,O] (8) 



for small k around the T point. The corresponding spin 
textures are shown schematically in the inset of Fig. 3a, c 
for p x orbitals in upper and lower Dirac cones and in 
the inset of Fig. 3b, d for p y orbitals. Unlike p z or- 
bitals, the amplitude of p x and p y orbitals for the sur- 
face states is not isotropic, but has 20^ angular depen- 
dence around the Fermi surface, as shown in Fig. 4c. 
We may take the difference of the amplitudes between 
p x and p y orbitals, as shown by colors in Fig. 4a, b. 
Here more red means more p x character, and more blue 
means more p y character. The angular dependence in- 
dicates a tangential orbital texture for the upper Dirac 
cone and a radial orbital texture for the lower Dirac 
cone, as schematically shown by the inset of Fig. 4a, b, 
respectively. This orbital texture was experimentally ob- 
served recently 14 . Furthermore, we also plot the total 
spin textures for in-plane orbitals on the same figure, 
which show a right-handed texture for both upper and 
lower Dirac cones. All these salient feature can be un- 
derstood by the wavefunctions (4) and (5). For the up- 
per Dirac cone, although both | te)\p r ) and | le)\pt) 
terms exist in the wavefunction (4), their associated coef- 
ficients are unequal. When ^2 a (va. a — u\_ a k + w\_ a k) 2 > 
X)q( u o.q — ui ta k — u>i. Q fc) 2 , | le)\pt) term dominates over 
I te)\Pr) term, dominantly giving a tangential orbital tex- 
ture coupled to a right-handed spin texture. Similarly, for 
the lower Dirac cone, when ^2 a (vo, a + ui^ a k + wi^ a k) 2 > 
J2a( v o,a + ui^ a k - wx^k) 2 , \ U}\Pr) term in the wave- 
function (5) is dominant, yielding a radial orbital texture 
coupled to a right-handed spin texture. The difference 
between p x and p v orbitals can be calculated directly 
as D p*,o - D p v ,o = T2 cos 26> fc [(v , a T ku lta )kw lja ], 
which indeed shows a 29k angular dependence, and 
the total spin textures for in-plane orbitals can be ob- 
taincd as [£>±,I>±] = [£>± , x + D± x ,D± y + D± J = 



4[— sinflfc, cos#fc] J2 a ( v o,a T kui^ a )kwi, a , which shows a 
right-handed spin texture when ^2 a {vo ya ^kui ya )kwi^ a > 
0. Especially, if k gets close zero, both the total spin tex- 
ture [D^,Dy] of the in-plane orbitals and the difference 
between andp y orbitals — D^ approaches zero, 
also as shown in Fig. 4a, b. 

Therefore, there is a transition from a tangential or- 
bital texture in the upper Dirac cone to a radial orbital 
texture in the lower Dirac cone, switching exactly at the 
Dirac point. To quantitatively describe this transition, 
we introduce a polarization quantity 



p p A±) = 



£> Pj ,o(±, = 0) -£^,(,(±,0 = 90) 
D p .,o(±,0 = O) + D p .,o(±,0 = 9O) 



(9) 



with '±' for upper and lower Dirac cones. The plot of 
P Px (±) is shown in Fig. 4d where the energy level of the 
Dirac point is shifted to zero. In order to show the plot 
more visually, we reverse the value of the P Px (+) for the 
upper Dirac cone plotted with the red. The feature of 
P Pa; (±) undoubtedly indicates that the state of the lower 
Dirac cone forms a radial orbital texture, and the state of 
the upper Dirac cone forms a tangential orbital texture. 
The Dirac point is shown to be the exact transition point 
from the tangential to radial orbital texture. This is ex- 
actly the behavior observed in a recent experiment and 
explained within the first principle calculations 14 . The 
numerical results fit well to the analytical calculation, 
with the expression 

2E tt (»o ; « T Eu 1 . a /hv f )Ew 1 , a /hv f 



p p A±) = t- 

E« [(«b,a T Eu ha /hv f ) 2 + E 2 w 2 Jh 2 v 2 

(10) 

with the energy E. For small E around Dirac point, 
P Px (±) cx =p ^°™)2° shows the linear dependence 
on energy, as found in Fig. 4d. 

Although in-plane orbitals show different spin textures 
compared to p z orbitals, we stress the p z orbitals (50%) 
dominate the states near the Dirac point with p x and 
p y only around 30%. Therefore, the spin texture for the 
whole states show left-handed for the upper Dirac cone 
and right-handed for the lower Dirac cone, the same as 
that of p z orbitals, as well as the total angular momentum 
texture. 



IV. DISCUSSION 

In order to detect the spin texture of electrons, the 
spin-resolved ARPES technology has been developed 
by taking advantage of spin-dependent scatting pro- 
cesses and precisely measuring the magnitude of the 
asymmetry in the spin-dependent intensity with per- 
fect spin-polarimeters. 12 The non-trivial spin texture 
of surface states of TIs has been clearly observed by 
experiments. 8-12 In addition, the orbital character can 
be detected through the photon polarization selection 
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rules 22 based on the symmetry analysis. With this tech- 
nology, the orbital texture of surface states of Bi 2 Se3 was 
reported recently by a polarized ARPES experiment. 14 
Therefore, it is possible to combine these two technologies 
together in an electron spin-resolved and photon polar- 
ized ARPES experiment, with both the spin and orbital 
textures extracted in the same measurement. The pre- 
dicted spin-orbital texture can be directly confirmed in 
this type of experiment, which can explicitly reveal how 
SOC plays a role in the real material at the atomic level. 
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